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Abstract 

A large deviation principle is derived for stochastic partial differ- 
ential equations with slow-fast components. The result shows that the 
rate function is exactly that of the averaged equation plus the fluctu- 
ating deviation which is a stochastic partial differential equation with 
small Gaussian perturbation. This also confirms the effectiveness of 
the approximation of the averaged equation plus the fluctuating devi- 
ation to the slow-fast stochastic partial differential equations. 

1 Introduction 

Uncertainties (noise) is widely recognized in modeling, analyzing, simulat- 
ing and predicting complex phenomena |2j [18j E2 |35l e.g.]. Noise causes 
rare events in nonlinear stochastic system describing the metastability of the 
system (THl E21 ED]- The classic example is a tunnelling event between two 
stable points of a macroscopic system. This rare event eventually occurs 
after a long time scale by an addition of a small external noise to the macro- 
scopic system, but the probability of such rare event converges to zero as the 
strength of noise tends to zero. We need to understand the rate of such con- 
vergence. The theory of large deviations provide a powerful tool to give an 
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estimate to the rate of such convergence, shown to be exponential for finite 
dimensional stochastic systems [SHI 1201 EH e.g.]. To study the metastability 
of a macroscopic system with small noise we must build the large deviation 
principle (ldp). 

Stochastic partial differential equations (spdes) are appropriate mathe- 
matical models for many multiscale systems with uncertain influences [38J. 
The LDP for spdes has been studied by many people [71 [TT | [T3"l [T^ \W \ I2TI 
EH EJJ [321 e -g-] under a different framework to that used here. However, 
there are very few results on the LDP for nonlinear stochastic system with 
two widely separated timescales, as often appears in a complex system. Frei- 
dlin and Wentzell [20] first studied the LDP for finite dimensional slow-fast 
stochastic system with partial coupling. They used a bounded assumption 
on the nonlinearity. Then Veretennikov [37] built a LDP for the full coupled 
case with bounded assumptions on nonlinearity. Later still, under a bounded 
assumption, Ioffe studied the LDP for stochastic reaction-diffusion equation 
with rapidly oscillating random noise in the special case when there is no cou- 
pling between the slow component and fast component [2"3"l [21] • However, 
there appear to be no other previous ldps for slow-fast coupled stochastic 
partial differential equations. 

This article establishes the Freidlin- Wentzell LDP for a class of spdes with 
stochastic fast component and deterministic slow component. Let D be an 
open bounded interval and L 2 (D) be the Lebesgue space of square integrable 
real valued functions on D. Consider the following pair of stochastically 
forced, coupled, reaction-diffusion spdes for any e > 

d t u e = d xx u* + f(u% v e ) , u e {0) = u e L\D) (1) 




v c (0) = v Q G L\D) (2) 



with zero Dirichlet boundary condition on dD. Here W(t) is an L 2 (D) valued 
Wiener process defined on a complete probability space (Q, J 7 , P) as detailed 
in the next section. In physical applications this supposition is that the 
noise directly drives microscopic modes v e — the noise only emerges in the 
macroscopic modes u e through nonlinear coupling. 

If parameter e > is small so that v 6 fluctuates rapidly, then an effective 
approximated system is desirable. Under some appropriate assumptions one 
averages the slow part u e over the fast part v e which yields the following 
so-called averaged equation describing the dynamics of the system on a slow 
time scale 

d t u = d xx u + f(u) , u(0) = u , u\ dD = 0. (3) 
Here f(u) is the average of f(u,v) over the distribution of the fast part v. 
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Cerrai et al. [El El e.g.] recently developed more on the averaging principle of 
SPDEs. Wang and Roberts [5H] very recently gave a further approximation 
result via a martingale discussion which shows that the deviation u € (t) —u(t) 
is approximated by y/ez(t), for some Guassian process z(t), in the case that 
fast component is coupled with slow component and without any Lipschitz 
assumption on the slow component. Cerria [TU] obtained the same result for 
the special case where the nonlinearity is Lipschitz and there is no coupling 
of the slow component to the fast component. Then the deviation estimate 
shows that u e approximates u with a small Gaussian perturbation and this 
suggests a LDP for {u e } e . By studying the LDP for some auxiliary systems, we 
prove the LDP for {u e } e , Theorem |9j Moreover, the rate function for the LDP 
of {u e } e in the main result is exactly that of {u e } e solving (HO p - flTT]) . which 
is the averaged equation fl3j) plus deviation up to errors of 0(e), Section [2J 
Our results further confirms the effectiveness of the averaged equation plus 
deviation to approximate slow-fast SPDEs ([I])-®. 

Recently, a weak convergence approach, which avoids giving some techni- 
cal exponential tightness estimates, was applied to obtain LDP for SPDEs [T71 
UU e.g.]. But this approach does not work here because a drift transforma- 
tion leads the fast system to become a non-autonomous system for which one 
cannot average the slow part over the fast part. For this here we still give 
some exponential tight estimates, Section [31 and then by some contraction 
principles and an approximation we obtain LDP for some auxiliary systems, 
Section HJ An approximation shows that the slow-fast stochastic system (Trj)- 
d2J) is comparable with the auxiliary systems near some functions, Section [51 
and we derive the LDP for {u e } € . 

Section |6] presents an example slow-fast reaction-diffusion spde to illus- 
trate the LDP theory. Section [7] then explores the parameter regime near a 
stochastic pitchfork bifurcation in this example. Constructing the stochastic 
'superslow' manifold and the evolution thereon confirms that there is indeed 
a close correspondence between the original example system and the LDP 
averaged system. 

2 Preliminaries 

Let H = L 2 (D) with L 2 -norm denoted by || • ||o and inner product by (•, •). 
Define the linear operator A = d xx with zero Dirichlet boundary condition 
on D. Then operator A is the generator of a compact analytic semigroup e At , 
t > , on H. Moreover, denote by {ej}^ 1 , which forms a complete standard 
orthogonal basis in H, a family of eigenf unctions of A and — Aei = A^j , Aj > 
0,i = l,2,.... For any a > and u & H define \\u\\ a = \\A aj/2 u\\o ■ Then let 
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Hq be the space of the closure of Cq°(D), the space of smooth functions with 
compact support on D, under the norm || ■ || a . Furthermore, let H~ a denote 
the dual space of Hq . Also we are given H valued Wiener processes W(t), 
t > 0, defined on the complete probability space (ft, J 7 , Tt, F) [27]. Denote 
by E the expectation operator with respect to P. We consider the spdes of 
the form (JXJ) — (I2J) with separated time scale and with a ^ is an arbitrary 
real number parametrising the strength of the noise. We adopt the following 
four hypotheses. 

Hi f(x,y) : RxR — > R is Lipschitz continuous in both x and y with Lipschtiz 
constant Lf and for all x, y £ R 

\f{x, y) \ 2 < ax 2 + by 2 + c , f(x, y)x < ax 2 + bxy + c , 

for some positive constants a, b and c. 

H2 g{x,y) '■ M x R — )■ R is Lipschitz continuous in both variables with 
Lipschitz constant L g . For any x,y e R 

g(x,y)y < -dy 2 + exy 

for some positive constants d and e. 

H 3 L 9 < Ai . 

H 4 W is a Q- Wiener processes which has the following series expansion 

00 

W(t,x) = J2VQiei(xMt) ( 4 ) 

i=i 

and 

Q&i Qi&i 7 ^ 1,2,.... 

Moreover, tr Q < 00 . 

By the above assumptions we have the following results on the averaging 
approximation to the slow-fast spde ([I])-© [39|. 

Theorem 1. Assume i?2- For any fixed slow part u G H , the fast system (TJ|) 
has a unique stationary solution, r] €,u (t), with distribution fi u independent 
of e. Moreover, the stationary measure fi u is exponentially mixing. Also, 
rj e ' u is differentiable with respect to u with Frechet derivative 

D uV € ' u < D v 

for some positive constant D v which is independent of e and u . 
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Given the stationary measure /i 11 for the fast part, we define the following 
deterministic averaged equation 

du = [Au + /(«)] dt , (5) 
w(0) = u , (6) 

where the average 

/(«)= / f(u,v)n u (dv). (7) 
Denote by por the metric on space C (0, T; if) with 

Pot(u, v) = m ax^ ||u(t) — v(t)||o for all w , t> G C(0, T; ii) , 

then we have the following theorem. 

Theorem 2. Assume H\-H±. Given some T > 0, for any uq G H , 

solutions u e (t,u ) of (TJP converges in probability to u in C(0,T; H) which 
solves (EP-(G|). Moreover, the convergence rate is 1/2; that is, for any k > 

P{p 0T (w» <C^v^} >l-« 
/or some positive constant C£ > . 

Now define the deviation between solutions u e and averaged solution u, 

z ^ t )= * (u--u). (8) 
For e = 1 we write r? 1 '" = 77" . Then we have 

Theorem 3. The deviation z e converges in distribution to a stochastic pro- 
cess z in the space C(0, T; H) which solves the spde 



z = Az + f' u {u)z + y/B(ujW (9) 
where B(u) : H H is Hilbert-Schmidt with 



00 



B(u) =2 E [(/(«, ^(t)) - /(«)) ® (/(«, rf (0)) - /(«))] (it 
#(«) = / 



and is an H -valued cylindrical Wiener process with covariance opera- 

tor Id-H- 
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Remark 4. Assume {ei} c *L 1 be a standard eigenbasis of H, then B(u) has 
the following series form 



^ />oo 

%) = 2V / E </( M ,r ? u (t))-/»,e,)</( M ,r / u (0))-/»,e :; > rfte, <g> ej 
i,j=i Jo 1 J 

where e^ ® ej is the tensor product of e, and e^ . 

Then formally we write the following averaged equation plus deviation 
up to errors of (9(e) as 



du e = [Au e + fiu')] dt + y/ly/B{u e ) dW{t) , (10) 
M e (0) = u . (11) 

A family of random process {u € } t in space C(0,T; H) is said to satisfy 
the LDP with rate function / if 120 



1. (lower bound) for any ip G C(0, T; if) and any S, 7 > , there is an e > 
such that for < e < e 

V{p«r(u e ,<p)<8}>exp* /( ^" ) + " 

2. (upper bound) for any r > and any 5, 7 > , there is a e > such 
that for any < e < e 



F {p 0T (u% K T (r)) >S} <exp 



r — 7 



Here the level set lf r (r) := {y5 G C(0,T;H) : /(</?) < r} . If K T (r) is 
compact, then rate function I is called a good one. 

There are many known results on the LDP for those SPDEs in the form 
of (TXOT) — fTTTT) as e — > , as mentioned in Section [TJ Thus one may expect 
to derive LDP for {u e } e from ffTOl - ffTTl) . However, it is difficult to obtain 
an exponential approximation in probability between u e and u e as e — > , 
which is needed to pass the LDP of {u e } e to {u e } t . But here our result shows 
that the rate function for {u e } t as e — > is exactly that of {u e } e - This 
also shows system (fTU|) - ([TT|) is indeed an effective approximate description 
of the macroscopic behaviour of the slow-fast system ([I])-©, even when one 
considers the exit problem caused by small noise perturbation for u e in the 
system of slow-fast SPDE CO)-©. 

For our purposes we study the LDP for a series of auxiliary slow-fast 
stochastic systems by a contraction principle and an approximation argu- 
ment. Then a controlled approximation derives the LDP for {u e } e . We first 
give some contraction principles [15] which are used in our approach. 
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Lemma 5. Let X and y be Hausdorff topological spaces and <3> : X — > y a 
continuous function. Consider a good rate function I : X — > [0, oo] . Then 



1. for each y G y , 



J(y) = inf {/(*): 1/ = <&(*)} 



is a good rate function on y , where the infimum over an empty set is 
taken as oo; 

2. if I controls the LDP associated with a family of probability measures {fJ, £ } e 
on X, then I controls the LDP associated with the family of probability 
measures {ji e o on y . 

To introduce a generalized contraction principle we give the following 
definitions. 

Definition 6. Let (y, d) be a metric space. Then the probability measures 
and {ft e }e on y are called exponentially equivalent if there exists a prob- 
ability space {Q, £> e , P e } and two families of y~ valued random variables {Z e } € 
and {Z t } t with joint laws and marginals {ji e } and {ji e } e respectively, 

such that for each 5 > , the set {cu : (Z € , Z e ) £ Ts} is B e measurable, and 

limsup elogPefT^) = — oo 

where T s = {(y, y) : d(y, y) > 5} C y x y . 

As far as the LDP is concerned, exponentially equivalent measures are 
indistinguishable. 

Lemma 7. // an LDP with good rate function I(-) holds for the probability 
measures {fi e }, which are exponentially equivalent to {fi e }, then the same 
LDP holds for {fie}- 

Furthermore, we need the following assumption 

H 5 There is a positive constant c and c 1 such that 

(B(ip)h,h) > c \\h\\ 2 and (DB(ip)h, h) < Ci||/i||g , for all ip , h e H 

where DB is the Frechlet derivative of B . 

Remark 8. Under the above assumption, a/ 'B(<p) is Lipschitz continuous 
in ip. The following is one simple example of / such that B(ip) satisfies 
assumption H 5 , 

f(u,v) = f 1 (u) + f 2 (v) 
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with fx and f 2 both Lipschitz continuous. If the stationary Gaussian pro- 
cess rj u = 7] is independent of u, by the expression of B(u) in Theorem El 



B(u) = 2 E 



oo 







(f 2 ( V (t)) - Ef 2 (r,)) ® (/ 2 (ry(0)) - E/ a fa)) 



Then B(u), which is independent of u, satisfies assumption H 5 . Section [6] 
details one example where /2(f) is linear in the fast v. 

Now for the slow-fast stochastic system (JTJ we define the following skele- 
ton equation 

<p = A<p + f(<p) + y/B(p)h, <p(0)=uo. (12) 

And define the rate function 

1 ' T 



I u (<p) = inf ^ / IIZiOOHjjds : y» = </ (13) 

where ip h solves f )12p and with inf = +00 . Then we give our main result. 

Theorem 9. Assume H\ H§ . For any T > , {u e } e satisfies the LDP with 
good rate function I u on space C (0, T; H) . 

Remark 10. By the large deviations principle for stochastic evolutionary 
equations [27], under assumptions H4-H5 , {u e } e satisfies the LDP with rate 
function I u (ip) in space C(0, T; H) for any T > . So the averaged equation 
plus deviation (fT0]) - (fTT]) predicts the metastability of ([I])-©. 



3 Exponential tightness 

We prove some exponential tightness results which are crucial for the LDP 
estimates. 

For any ip € C(0, T; H) consider the following system 

du^ = [Au e ^ + f{ip, v^)] dt , u e '^(0) = u , (14) 

dv^ = - [Av^ + g(ip, v e ^)] dt + -^dWit) , u e '^(0)=-ub. (15) 
e V e 

Then we have the following exponential tightness result. 

Lemma 11. Fix ip e C(0, T; if). For any T > and e > , {u t '^}o<t<e is 
exponentially tight in space C(0,T; if). 
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Proof. This result follows via a uniform estimate to u e ^ in space C(0, T; i?o)H 
C a (0, T; #) for some a > . By ([II]) 



u e>v, (t)=«o+ / Au e ^(s)ds + / f{iP{s),v e ^(s))ds. 



Then by the increasing property of / 



\\u^{t)-u^{r)\\l 
< 2 17 ||At2^(s)|| tfo 



+ 2 



||/(V>(s),^(s))|| ^ 



< 2|f-r| / ||^( a )||gd S + 2|t-r| / ||/(^ S ),^(s))llo^ 



< 2|t — t| 
and 



T 



T 



\\Au^{s)\\lds + 2|f-r| / [a\\^(s)\\l+\\v^(s)\\l + c]ds. 



2dt n 111 



= -\\Au^\\l-(f^,v^),Au^) 
< -\\\Au^l + \\\f^^)\\l 



1 



1 



< -oll^HS + 5 [a|Hl2 + ||iJ e '1lS + c]. 



That is, 



||A^(5)||^5+ SUp 

0<s<T 



< 



l + a I U(s)\\'ods+ I \\v^(s)\\lds + cT 



Then to give a uniform estimate in space C(0, T; ii^) n C a (0, T; if) for some 
a > , we just need to estimate J \\v e ^ (s)\\q ds . 

Next we estimate J Q T Ht^Hi o?5 ■ The Lipschitz property of g and applying 
the Ito formula to H^'^Ho yield for some positive constants d and c" 

1 d 



~e,-0 II 2 



2dr 



\v 
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--ll^ll? + -(g&,v^),v^) + ^(^, wo + f trg 

< -jll^ll? + 7 IHIo + ^ tr Q + -L<^, WO • 
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Then 



e\\v 

ft 



< 2c" 



' < mi 2 + 2c' f \\v e 'f(s)\\lds -e\\v \\l 
Jo 

s)\\ 2 ds + tTQt + 2y/i f (v^(s),dW(s)) 

Jo 



Now define Ml = f*{v^{s),dW{s)) and A = c'A 1 /2g max with g max = 
maxj qi . And denote by {M)\ the covariance of M t e . Then we have 

c' J*\\v^(s)\\lds < ||,; || 2 + 2c // ^ < ||^)||^ S + trgt + v ^A^-^(M)^ 

+ ^-(M)t-c' J\\v^(s)\\lds 

2c" jf U(a)\\lda + tr Qt + V~eMt - ^<M)« . 



< KIIS + 

And by the exponential martingale inequality we have 

P l^eMt - ±{M)l > 5} = P \^MI - f e {M)\ > ^} < e-^ 

which yields the exponential tightness of {u e, ^} e . □ 

Similar method yields the following result. 
Lemma 12. For any T > , {u e } t is exponentially tight in space C(0, T; H). 

Proof. Follow the proof of Lemma [TT| we just need to estimate J Q T 1 1 § cZs-(- 

lo \\ vt { s ) 1 1 d s ■ By our assumption we estimate 

\\u e {s)\\\ds+ [ \\v e (s)\\ids. 
Jo 

For this, applying the Ito formula to ||f e ||o and by the assumption on g we 
have 

Moreover, 

\l t \\u%<-\\u%\-a\\u%-rb{u^y)+c. (17) 
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T r T 
2 



Then + {TTD yields 

|w e (s)||^s + 2 / \\v e {s)\\{ds 
Jo 

< \\uo\\ 2 +\\v \\ 2 + C T + trQT + y/~e I (v e (s) , dW (s)) . 

Jo 

Then a similar argument as for Lemma [TT] yields the exponential tightness. 
The proof is complete. □ 

4 LDP for some auxiliary systems 

Next we study the LDP for some auxiliary systems. 

First, let ip E H , which is independent of time, and consider the following 
system 

du e ^ = [Au e >^ + f(?P,v e ^)] dt, u € ^(0)=u , (18) 

dv^ = - \Av^ + g(if>, v e >^)) dt + —dWit) , v e '^(0) = v . (19) 
e V e 

We study the LDP for {u e, ^} e . For this define 

P+{t)= f [f{^y^{s))-m] ds. 



Then 

if>* = Au e '* + f(iff) + , u e '^(0)=u . (20) 

Now in order to obtain the LDP for {u e, ^} e we study the LDP for {u e, ^} e which 
solves 

fr* = Au^ + f{$) + i £ ' f , u^(0)=u (21) 

with 



Jo 



where rf'^ is defined in Theorem Q]. Now define the skeleton equation to (J2T 
with ip £ H 



= Aip f,h + + ^ ; ^(o) = Mq (22 ) 
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for h G L 2 (0,T; H) . Equation (|22|) is a linear equation, so for any ij) G H 
and /i G L 2 (0,T;H) there is a unique ^ G C(0,T;H) solving Q22]). Then 
define the functional 



We prove the following LDP result for {u e, ^} t . 

Theorem 13. Fix ip G H . For any T > , {u e '^} e satisfies LDP in space 
C (0, T; H) with a good rate function 1$ . Moreover 1% is lower semicontinuous 
in ip G H . 

Proof. First by the same discussion in the proof to Lemma [TT], for any T > , 
{u £ '^} e is exponential tight in space C(0,T; H). Then there exists {Kr\r 
which is a nondecreasing family of compact sets such that 

¥{u e ^ eK R }>l- e~ R/e . (23) 



Now for any fixed N > 1 , let if at := span{ej : i = 1,...,N} and P, 



N 



H — > Hn be the orthogonal projection. Then := PnC'^ e C(0,T; H 



N) 



satisfies LDP with good rate function [20 



iIn(Vn) = inf |i \\h N {t)\\ 2 HN dt : ^/B N (ip)h N = (p N ^ 
where ip^ G C(0,T; H N ) and 

POO 

B N {^) = 2 / E [(JWM ^ (*)) " *W(V0) ® (*WM ^(0)) - PnM)] dt . 



Introduce process uffl solving 

Ujf = Aujf + f N (ip) + t^f , «1?(0) = PnUq 

with /jv = Pjsrf . By the continuity of the map £jy H- u e jf in space C(0, T; H) 
and the contraction principle Lemma[5], {u e ^} t satisfies LDP with a good rate 
function 

itM = heL ^ T;H) \\Ks)\\lds : <p N = # 

where ip^' h G C(0,T; if) solves the following equation 

<p%* = A<p% h + f N (i(>) + y/BM$h, ^> h (0) = P N u 
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Moreover ip^ h — > tp^ ,h as N — > oo in space C(0, T; if) . Then we have 



I+(<p) = lim /J iJV (Piv^) = sup ii^Pstp) , <p G C(0, T; if) (24) 

which is a good rate function. 

Now for any 7 , 5 > , cp G C(0, T; H) and ff > , for u £ '^ G K R , there 
is N R (6, <p) > such that 

Por(u e ^, + Pot(<P, Pjv¥>) <S/2, N> N r (6, <p) . 

Then there is eo > such that for < e < eo 

F{p 0T (u^, V )<5} (25) 



> ¥\poT(uf,P N ip)< 5/2, N>N R (5 } <p) 



G K 



R 



> 



> 



> 



P [p 0T (u e /, Pn^p) < 5/2} P {u^ G K r I P ot(^, P^<p) < 5/2 

ii N {p NV )+i 



cxp 



exp 



P { G 



Poriujf , P N (p) < 5/2 



tf(¥>)+7 



P <! G if 



Pw{v$,P N ip)<8/2 \ . 



Passing to the limit R — > 00 and noticing that lim^oo P{w e '^ e i\^} = 1 , 
the above estimate yields the lower bound estimate. 

Now for any r > and 7 > , there is an e\ > and R > such that for 
any < e < t\ 



Pjw^ G < exp 



r — 7 



(26) 



Let Kt(t) = {y? G C(0, T; if) : f^f (<p) < r} and # TiiV (r) = G C(0, T; 

^wJv(VJv) — r ) • Then for w e '^ G if h , there is N(R,r) < such that 
for iV > N(R, r) 

¥{p 0T (u^,K T (r))>5\u^ eK R } 

< P {pot(«^, #iyv(r)) > 5/2 | ^ G 
r — 7 



< exp | _1L (27) 

for e is small enough. Then ( I2"6"|) - (j2"7j) yield the upper bound estimate. 
Next we prove the second result. This is followed by proving the set 



*(r) := {^eH:lt^)<r} 
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is closed for any r > and (p G C(0,T;H). For this proof, let {4> n } n C 
\l/(r) and ^ n — )■ ip in space if. Moreover, there is a sequence /i„ G £ 2 (0, T; if) 
such that (p^ n ' hn = (p with 

<p* n ' hn = Aip^ h - + /(^ n ) + y r W^n)h n , yA'' 1 " (0) = u 

and 

i / \\h n (s)\\ 2 ds<r+^. 
z Jo n 

Then there is a subsequence, which we still denote by h n , weakly convergent 
to some h G L 2 (0, T; if) . We show that (p^' h = up with 

0<M = A(p *,h + ^ + /b^)/, , ^,h( ) = Mo 

and 




By the assumptions Hi and H 5 , we have 

/ty„)->/M in H 

and 

VWn) v / W) in £(#). 

Then 

\/B{i> n )h n -> v^)/i weakly in L 2 (0, T; if) . 
The proof is complete. 

□ 

By the above result we have the following corollary. 

Corollary 14. For fixed ip G H , {u e '^} e satisfies LDP with good rate function 
1$ in space C(0, T; if) /or any T > . 

Proof. We show that -u e '^ is exponentially equivalent in probability to -u e '^ . 
Let £7^ = u e ^ - u e ^ , then 

U e >* = AU e '^ + t'^ - t'^ , U^(0) = 0. 

Notice that 

?>*(t)- ?*(t)= f [f{i>y^{s))-f{i>^{s))] ds. 

Jo 
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Then for any T > 
sup \\u^{t)-u^{t)\\l < C Xl sup f\\f^y^{s))-f^^{s))\\ Q ds 

0<t<T 0<t<T Jo 

< L f C M f \\v^(s)- V ^(s)\\ ds (28) 
Jo 

for some positive constant C\ 1 which depends on Ai . Let ( e '^ = v e ^ — rf^ , 
= -At?* + i [g(^ v e >*) - rj^)] , C'HO) =v - 7^(0) . 

Then 

By assumption H3 we have 

IIC^(*)ll^<e- 2te / e |ko-^(0)||^ 

for some positive constant c which is independent e . Then by fl28|) we have 
for any T > 

sup \\u^(t)-u^(t)\\ < L f C Xl [ T e- s ^\\v - V ^(0)\\ ds 

0<t<T Jo 

< etflh-^Uo 

for some positive constant C which is independent of e . 
Then we have for any 5 > , 

p{ sup \\u e ^{t)-u e ^(t)\\ > 5} <P{eC||w -?7 e,V, (0)|| > 5} . 

{0<t<T J 

And by the Guassian property of t] e, ^(0) , we have for any 5 > 
limsupelnpj sup \\u £ ^{t) - u 6 ^(*)|| > s\ = -00 . 
Then the generalized contraction principle Lemma [7] completes the proof. 

□ 

Now we consider the special case that ip is a step function which we denote 
by ip n . Let (u e * n , v e * n ) solve fn4|) - fTT5|) with ip replaced by tp n , we show that 
{u e * n } e satisfies the LDP with a good rate function. Typically we choose 

n-l 

1> n (t) = Y,WXluM (29) 

i=0 

15 



with i/)? G H and ti = iT/n , i — 0, 1, . . . , n — 1 . Let u e ^ and ■D 6 '*'* be the 
restriction of u e '^ n and v €, ^ n on the time interval £j + i], respectively, then 
we have 

which satisfies the LDP in space C([£j, t i+ i], H) with rate function I u i by 
Corollary [Ml Then we have the following result 

Theorem 15. For any n and step function ip n G C(0, T; if) in i/ie /orm / flP|) . 
{■u e '^ n } satisfies LDP iwi£/i good raie function 

where uo^ n,h solves UZ^) with ip replaced by step function ip n . 

Proof. For any up G C(0,T;H) with if (<p) < oo , let ^(t) = ¥>(t)X[ ti , tm] G 
C([ti, ti+i]; i?]) , then for any 8, 7 > 0, there is an e > such that for any 
< e < e , 



P <J max \\u e '^(t) - <pi(t)\\ < 8 } > exp 

ti<i<ii_l_! 



Moreover, by the assumption on /, for n large enough for any 8 > 

max \\u e >W(t) - Mt)\\o <s)d { \\u e '^(ti) - MU) Mo < 5' (8)} := A t 
for some 8' (8) small enough. Now for fixed time ti, u e, ^"(ti) is an if- valued 



16 



random variable. Then for small enough 5' = S'(S) 
P{pot(^»<5} 

max \\u e <W(t)-<pi(t)\\o<5 

l<i<n, ti<t<ti+i 



> P ( max \\u e ^(t i )-tp i (t i )\\ <5'\=F{A 1 A 2 ...A n } 

I l<i<n J 

= P {A x } F{A 2 \A 1 }¥ {A 3 | ^ A 2 } ■■■¥{A n \ A X A 2 . . . A n ^} 



> P <{ max ||u e >^°(t) - (p (t)\\o < 5' 



xP| max \\u e 'W(t) -<p x (t)\\ < 8' \ A x 



X ■ • • X P 



max 

in— 1 < ~~-i < ~~-i'n 



-i{t)-ip n - 1 (t)\\ Q <5'\A x A 2 ...A 



n-l 



i=l 



exp 



Now for any r > and 5 > , 

¥{p 0T (u^\K T (r))>6} 
< p{pot(« £ '^, # Mi+l] (r)) > 5 for some < i < n - l} 
r — 7 



< 



exp 



The proof is complete. □ 

The lower semicontinuous property of 1$ for any ip e C(0, T; if) is needed 
in our approach. We have 

Lemma 16. For any T > 7 I^if) is lower semicontinuous in both(p andip G 
C(0,T;H). 

Proof. The lower semicontinuous property in ip is followed by a similar dis- 
cussion in the proof of second part of Theorem [13] by the fact that both f(ip) 
and -\fBijp) are continuous in tp in space C (0, T; H) . And the lower semi- 
continuous property in ip is followed by the same discussion for the usual 
evolutionary equation [IT] . □ 

For every E C C(0,T; H), denote by Int(-E) the interior of E and Cl(E) 
the closure of E. Then by the above result and the lower semicontinuity 
of in if we then have, for any step function ip 6 C(0,T; H) in the 

form (EHD- 
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Corollary 17. For any E C C(0, T; if) , 

- inf rf(» < liminf elnP{u e '^ G -E} 

¥>eInt(B) e->0 

< limsupelnP-{>/^ G £} < - inf . 

Now for our purpose we need the following result 

Theorem 18. For any ft G L 2 (0,T; H) , the skeleton equation [Wjl has a 
unique solution ip G C(0, T;H). 

Proof. By the assumptions on / both f(<p) and y/B(ip) are Lipschitz con- 
tinuous. Then the result follows by a standard discussion on deterministic 
pdes [TTJ . □ 



Then we define the rate function I u (ip) = I£(tp) for any ip G C(0, T; H). 
Furthermore, for any ip G C(0, T;H), the following relation between rate 
functions 1$ and I u is needed to derive the LDP for {u e } 6 . 

Lemma 19. Let ip n be a family of step functions uniformly converging to ip G 
C(0, T;H) as n — )■ oo . Tften tftere zs a family of functions ip n G C(0, T; if) 
converging to ip, such that 

limsup/f (¥>")<!„(¥>). 

n— >oo 

Proof. Suppose i u (v?) < r < oo . Otherwise the result is clear. 

By the definition of i u (<£>), there exists a function ft, and a sequence h n G 
L 2 (0,T; if) such that 



<p{t) = S{t)u + / 5(t - s)f{<p(s)) ds + / v / ^^(s))ft(s) ds (30) 
and h n -> ft in L 2 (0, T; if) with 

i / r ||ft"( s )||^ s </ u (^) + -. 

z Jo " 

Then define 

^ n (t) = S(iK,+ / S(t-s)f(i/> n (s))ds+ [ ^/B{^{s))h n {s) ds. 



By the Lipschitz property of /, the definition of \jB(ip) and that %p n — >■ 
in space C(0, T; if) we have 

y?" -> y? inC(0,T;ff) asn-)-oo 

and 

limsup/f M< 

n— >oo 

This completes the proof. □ 
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5 LDP for slow-fast stochastic partial differ- 
ential equations 

Now we show that {u e } e satisfies the LDP with a good rate function I u . To do 
this we prove a special relationship between u e and u e ^ in space C(0,T; H). 
The relationship shows that u 6 ^ is comparable with u € near tp and implies 
the LDP. Here we do not restrict ip G C(0, T; H) to be a step function. 
For this let = u e - u e >^, then 

U^ = AU^ + f{u e ,v e )-f{^,v e ^), U € '^{0) = 0. (31) 

In the mild sense 

U^it) = [ S(t - s) [/(««(*), v<(s)) - f^(s),v^(s))] ds 
Jo 

= f S(t - a) [f(u%s), v*(s)) - f&(s), v*(s))\ ds 
Jo 

+ [ S(t-s) [f^(s),v^s))-f^(s),v^(s))] ds. 
Jo 

By the assumptions on / and the analysis on fast motion v e ^ 

\\f{u*,v*)-f{Tp,v*)\\ < L^K-VHo, 
\\f(i/>,v e )-f(i>,v^)\\ < L f \\v £ -v £ '%<L f D v \\u e -i;\\ . 

Then there is a positive constant Ct such that 

\\u e - U^Wc^T-H) < C T \\U E - 1p\\c(0,T;H) ■ (32) 

Rewrite (l3Tj) as 

= AU e ^ + f{u e ^,v e ^)- f(tP,v e ^) + f{u e ,v e )- f{u e ^,v e ) 
+ f(u^,v e )-f(u^,v^) 

Then by the Lipschitz property of / we have one positive constant L depend- 
ing on Lf and D v such that 

^il^'llo < ^ll^'llo + L f (l + D v )\\u^ - m 
which yields that 

\\U^ - « e ||c(0,T^) < C T \\U^ - nC(0,T:H) (33) 

for some positive constant Ct- Now by relations fl32]) and f )33|) we prove the 
main result. 
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Proof of Theorem^ We follow Freidlin and Wentzell's [20] approach to ob- 
tain LDP for slow-fast random ordinary differential equations. For any (p £ 
C(0, T; H) with I u {<f) < oo and any 7, 5 > , by Lemma [191 we can choose 
a step function ip n and a function ip n such that 

p 0T (^»<i, max ||^(t)-^)|| <i and if (^) < I u {y) + I . 
Now by ([31]) 



< 


Pot( 






+ P0T« 




) +Pot(<p",v) 


< 


Pot( 






+ C*r Port 




<T) + Pot(<p", <p) 


< 


Port 






+ CtPotI 







+ C T max - ^(f) ||o + p or (^, <p) 

< (l + C T )p T(^,<p") + (2C T + l)i. 

n 

Then for any 5 > and 7 > , by Theorem [T5l there is an e\ > such that 
for any < e < ei we have the following lower bound estimate by choosing n 
large enough 

P{ P 0T«<P)<5} 

> F {(l + C T )p T(u^ n ,ip n ) <5/2} 

r /fv^+71 ( 4M + 27 

> exp < > > exp 



Now we prove the upper bound estimate. First, by the exponential tight- 
ness of {u e } e in space C(0,T; H), for any r > there is a compact set K r 
such that 

lim sup e In P {u € £ K°} < — r . 
Then for any 7 > there is an e 2 > such that for any < e < e 2 

r — 7 



P {u e £ i^} < exp 

As K r is compact, choose a finite <5'-net in K r with 5' < 5 and let (pi, (p2, ■ ■ ■ , <p r 
be the elements of this net, not belonging to Kt{t). Then 

n 

P {p 0T (u e , K T (r)) > 5} < F iPoT(u e , Vl ) < 5'} + P {u £ £ ^} 

i=l 
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Now we choose step functions ipi, ip 2 , ■ ■ ■ , VVi sucn that 

Parish <Pi) < , i = l,2, ...,n. 
Then by the inequality (13"31 for z = 1, 2, . . . , n 

¥{ PoT (u\ Vi ) < 5'} < p{p T(£ e *,<A-) < 2(C T + 1)5'} . 
And by Corollary [T7] for i = 1, 2, . . . , n we have 
p{por(fi e ^i) <2(C T + 1)5'} 
< exp {-- [inf : Po t(^, fft) < 2(C T + 1)5'} - 7 ] 

By the semicontinuity of the functional 1% in ip, Lemma [TBI, that for any 7 > 
, there is 5' such that 

It i (ip)>r -7/2 for pvr{<P,<Pi) <2(C T + 1)5' and Jf(^)>r. 

Then by the choice of ipt ^ Kx{r) we have 

p{por(« e ^S ifi) < 2(C T + 1)5'} < exp 
for e small enough. The proof is complete. □ 

6 An example of slow-fast stochastic reaction- 
diffusion equation 

Next we consider the following slow-fast spde on the domain (0, L) with zero 
Dirichlet boundary condition 

d t u e = d xx u e + A sinw e — v e , u e (0) = u , (34) 
e d t v e = d xx v e -v e + u e + y/ead t W{t) , v € (0) = v (35) 

where W is L 2 (0, L)-valued Q-Wiener process and A, a > are constants. As 
usual, the small parameter e measures the separation of time scales between 
the fast modes v and the slow modes u. For small e our LDP theory gives 
the spde (j37p as an appropriate (weak) model for the stochastic dynamics 
of the slow modes u. 

Apply the LDP theory to this system. Note that the nonlinear reaction/ 
interaction function f(u, v) = A sinu — v is Lipschitz continuous. Denote the 
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operator A = d xx with zero boundary condition on (0, L) . Now for fixed u 
the fast system (1551) has a unique stationary solution r) e ' u with distribution 

^ u =M[{I-A)-^a ^ I ~ A ^ Q 

Then 

f(u) — A sin u- (I - A)~ 1 u . 
Let 7] u be the stationary solution of 

d t v = d xx v -v + u + ad t W(t) 

for any fixed u G L 2 (0, L) . Then rf distributes as \i u and 

/■oo 

B{u) = 2E / - (J — v4) _1 m] ® [77^(0) - (/ - A)~ 1 u\dt 

Jo 

Noticing that 

E V u (t) ® (0) = a 2 exp {-(/ - A)t} 
we then have 



{I-A)-'Q 



+ (I — A) 2 u(g)u 



s/B$j) = {I - AY^ay/Q 

which is independent of u and satisfies assumption H 5 . By Theorem [HI 
{u 6 } satisfies LDP on space C(0,T; H) with good rate function 

T \\h{s)\\lds:v = V h 



= inf 

^' hdL\0,T;H) \2 

with ip h solving 

(p = A(p h + Asiny? - (J - A)~V + ( J ~ AY 1 cr^/Qh, tp(0) = u ■ 
Furthermore, the rate function is 



/-A 



[ip(s) - A^(s) - A sin ip + (/ - AY V] 



ds (36) 



for <p is absolute continuous. Otherwise I u {(p) = oo . 

Now we write out the averaged equation plus the deviation for flM|) - fl35]) 

as 

du e = [Au e + Asinu e - (J — A}' 1 ^] dt + yfea(I - AY 1 y / QdW(t) (37) 
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slow-fast spde 



superslow 
SDE flU 



dMD (ESD 



stochastic centre manifold 



LDP theory, 
sections [2HS] 




slow model 



superslow 
SDE flU]) 



spde (137j) 



stochastic centre manifold 



Figure 1: this schematic diagram shows how the construction of stochastic 
slow manifold models, in Section [7J verifies by the approximate equality of 
the two superslow models that the LDP theory correctly captures the slow 
dynamics of the slow-fast spde fl34|) - fl35|) . 

Then by the LDP for stochastic evolutionary equation [57] . {u € } satisfies LDP 
with rate function I u ((p) defined in (I36I) . This shows that the averaged equa- 
tion plus deviation (I3"7j) does describe the metastability of {u e } e solving flMty 
for small e. Moreover, for large enough parameter A, the spde model (I3"7j) has 
two stable states near zero for e = . When e ^ , noise causes orbits near 
one stable state to the position near the other one which shows the metasta- 
bility of the system fl37j) . So the slow-fast stochastic fl34|) -f l35|) also has such 
metastability described by system f l3"Tj) . The description of such tunnelling 
of the orbit needs detail analysis by the LDP which is left for future work. 

7 Stochastic centre manifold models confirm 
the LDP example 

This section uses the example slow-fast spde fl34|) - fl35|) to verify the LDP 
theory in a significant parameter regime. Without loss of generality set the 
non-dimensional domain length L = it. In the absence of noise the spde fl34|) - 
undergoes a deterministic pitchfork bifurcation from the trivial field u = 
v = 0, as parameter A crosses the critical value 3/2, to two nontrivial fixed 
points u ~ a/A — 3/2 sins . Consequently, with noise, a stochastic pitchfork 
bifurcation takes place in the vicinity of parameter A « 3/2 [3j |33j HJ e.g.]. 
As shown schematically in Figure HJ here we establish that the stochastic 
bifurcation dynamics of the original slow-fast spde ( I34p -( 131)|) and that of the 
LDP averaged spde ( 13"T|) are identical to the expected 0{e) error. 

As indicated in Figure [H we make a wide ranging comparison of the dy- 
namics by constructing and comparing the stochastic centre manifolds, and 
the evolution thereon, of the two spde systems near this stochastic bifurca- 
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tion. We explore dynamics in the vicinity of the stochastic bifurcation by 
setting the parameter A = | + A' for any small enough bifurcation parame- 
ter A'. For small A' and small e there are three time scales in the example 
fast-slow spde (|34|) -( 135|) : the v modes quasi-equilibrate on the fast time scale 
of 0(e); almost all of the u modes evolve on the slow time scale of 
but the sinx mode in u evolves on the superslow, long time scale of 0(l/\'). 
The LDP averaged spde has just the latter two time scales in this parameter 
regime. The interactions among three time scales is a major complicating fac- 
tor in constructing the stochastic centre manifold and the evolution thereon. 
Because of the three time scales, and to be consistent with the terminology 
of earlier sections, we henceforth refer to the stochastic superslow manifold 
as it is the evolution on the superslow, long time scales of C(l/A') that we 
encompass and compare in this section. 

The stochastic superslow manifold (ssm) is based from the linear dynam- 
ics exactly at critical and, for simplicity, based from no noise [29], CD e.g.]. 
Exactly at critical, and with no noise, A' = a = 0, both the SPDEs, (134"]) - 
([351) and spde (1571) . have centre subspaces about the origin: u = asinx, 
v = ^asmx and u = asmx, respectively. Stochastic centre manifold the- 
ory [TJ e.g.] then asserts that in a domain of small but finite amplitudes 
a and a, and small but finite noise a and small but finite parameter A', 
there exists an emergent SSM: solutions are attracted to the SSM roughly 
as exp(— |^t), and then evolve on the superslow long time scale. We compare 
the slow-fast spde (JH3D (JHSD with the LDP averaged spde (jHTJ) via construc- 
tion of their SSM models. 

Constructing ssms and the evolution thereon has many technical chal- 
lenges reported in detail elsewhere [121 ESI EH HOj e.g.]. The technicalities 
are even more challenging here due to the three times scales in the slow-fast 
SPDE (I34p - (l3l)j) when, as we assume, the parameter e is small. The construc- 
tion procedure used herein is detailed in a separate technical report [30J that 
all can check, reproduce and perhaps modify to other problems in the same 
class. For our purposes we appeal to a little more of the theory of stochastic 
centre manifolds: Arnold [I], building on the work of Boxler [6], assures us 
that if the SPDEs are satisfied to some order of residual in the small param- 
eters, then the ssms and the evolution thereon are constructed to the same 
order of error. Thus one may confirm the veracity of the following ssms by 
substituting the expressions into the SPDEs and verifying that the residuals 
asymptotically small as required. 

To reduce complicating detail but retain significant information in the 
example, we truncate the noise to its first three spatial modes: W = tpi sin x+ 
ip2 sin 2X+993 sin 3x where ipi denote formal derivatives of independent Wiener 
processes. Including more noise modes appears to just greatly increase detail, 
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without adding any significant change to the nature of the interactions seen 
among these three modes. 



Slow-fast SPDE flM|) — (j3"5"|) In six iterations, computer algebra [SU] con- 
structs the stochastic superslow manifold model. In terms of the superslow 
evolving amplitude a(t), where u ~ asinx and v ~ ^asinx, the stochastic 
bifurcation sde for the amplitude is 



16 + 8 A + 64 6 ) a 



a = A'(l + \e)a 

- v 7 ^ [(§ + y)<fi + i^Vs] 



27, 



^a 5 

9728 



.38 , 



+ ecr a 
+ C(a 6 + A /3 + e 3 + a 6 ) 



6080 



(3? 



In this and other expressions, convolutions e at * <p = J °° e a,s ^(t — s) ds . 
The corresponding SSM has slow field 



u = asinx + -^a 3 sin3x + ^y/ea sinx e e **y?i 



27, 



38 , 



±sin2x(e io c *-e e 1 * )ip 2 + ^ sin 3x( e 5 1 * - e s 1 *)^ 



+ 0(a 4 + A /2 + e 2 + ( x 4 ). 
and the fast field 



(39) 



v — I a sin a; + Y2T6 fl3 sm ^ 



O" 

H -= sm x 



2 22 
(1 + ±e)e~ 6* * +|e"e* * e~ 7** 



Pi 



(7 



H — p sin 2x 



5 27 5 5 



(7 



H — -p sin 3x 
+ 0(a 4 + A' 2 + e 2 + a 4 ). 



^2 

L _io 1 -¥*^ , 1 *, -±V 

Too e ) e e *- e Too e 5 *+Io e € * e e * 



^3 



(40) 



Zr 

This fast field v has large fluctuations, through terms like ~^ e ~ 6 be- 
cause convolution over the fast time scale, e - ^*, is O^^/e). However, note 
that repeated such convolution, e~ /3i//e *e~ /3t//e *, is C(e 3 / 2 ) [281 equation (27)]. 
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LDP averaged SPDE (|37|) In just four iterations, computer algebra [30] 
constructs the stochastic superslow evolution to be 



+ eo~ 2 a 



+ C(a 6 + A /3 + e 3 + a 6 ). (41) 
The corresponding SSM is 



m = a sin x + g§g« 3 sin 3x 



i -2^ t 1 

~ sin 2x e 10 * </> 2 + 777 sin 3x e 5 * ip 3 



- y/ecr 

+ 0(a 4 + A /2 + e 2 + a 4 ). (42) 

Compare the two superslow models First compare the slow field u for 
the slow-fast SPDE, ([39]), with the slow field for the LDP averaged SPDE, (j4"2"jh 
The differences are the 0(e) terms in the fast time convolutions ^/ee - ^*^*. 
Since u is u averaged over fast fluctuations, these differences are acceptable, 
and also ensure that the two amplitudes correspond: a — a + 0(e). 

Second, compare the evolution of the amplitudes, (1381) and (jHJ). The only 
differences are in terms 0(e), as indicated schematically on the right-hand 
side of Figure [TJ Thus the dynamics of the two superslow sde are within the 
claimed accuracy of the LDP. We conclude that this section verifies that in a 
parameter regime at least near the stochastic bifurcation, the LDP averaging 
approximation is correct. Although the LDP averaging only assures us that 
the slow model is correct in a weak sense, the strong identity between the 
convolutions appearing in the two SSMs, fl38l) - fl39|) and ([4"I]) - (j4"2"]) . suggests 
the correspondence between the LDP averaged system and the original system 
is generally stronger. 
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